Cu 3 Nb 2 08: a multiferroic with chiral coupling to the crystal structure 
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By combining bulk properties, neutron diffraction and non-resonant X-ray diffraction measure- 
ments, we demonstrate that the new multiferroic Cu3Nb20s becomes polar simultaneously with the 
appearance of generalised helicoidal magnetic ordering. The electrical polarization is oriented per- 
pendicularly to the common plane of rotation of the spins — an observation that cannot be reconciled 
with the "conventional" theory developed for cycloidal multiferroics. Our results are consistent with 
coupling between a macroscopic structural rotation, which is allowed in the paramagnetic group, 
and magnetically-induced structural chirality. 

PACS numbers: 75.85.+t, 61.05.F-, 61.05.cp 



Magnetic multiferroics, where the electrical polariza- 
tion emerges at a magnetic ordering temperature due 
to symmetry breaking by the magnetic structure, have 
been intensively studied in the past few years. The 
vast majority of these new materials, such as TbMnO.3 
P, i], Ni 3 V 2 8 Q CoCr 2 4 0, and MnW0 4 @ are 
cycloidal multiferroics, so called because their magnetic 
structures can be described as incommensurate circu- 
lar (or elliptical) modulations with the wavevectors in 
the plane of rotation of the spins. In cycloidal multi- 
ferroics, the non-collinear magnetic configuration itself 
is established by the competition between nearest and 
next-nearest neighbor interactions. The coupling to the 
crystal structure occurs through the spin-orbit interac- 
tion, making it energetically favorable to develop local 
Dzyaloshinskii-Moriya (DM) vectors, associated with a 
local polarization. In simple, high-symmetry cases, the 
electrical polarization is perpendicular both to the mag- 
netic propagation vector and to the normal to the plane 
of rotation of the spins, since the following formula holds 
Q: P — Xk m x (si x s 2 ), where Si and s 2 are two ad- 
jacent spins along the propagation direction k m and A is 
a coupling constant. In more complex, lower symmetry 
cases, the electrical polarization need not be perpendic- 
ular to the propagation vector, which, in turn, need not 
be contained in the plane of rotation of the spins (generic 
helicoidal structures). However, it is a strong prediction 
of the cycloidal multiferroics model that if all the spins 
rotate in a common plane, then the electrical polarization 
must be strictly contained within that plane. 

In this letter, we present a new multiferroic, with 
chemical formula Cu3Nb 2 08 and centrosymmetric tri- 
clinic symmetry (space group PI) in the paramagnetic 
phase. Using magnetic neutron powder diffraction, mag- 
netic susceptibility, heat capacity, electrical polariza- 
tion and non-resonant X-ray magnetic scattering mea- 
surements, we show that Cu3Nb 2 08 orders magneti- 
cally at Tn ~ 26 K, and develops an electrical polar- 



ization below a second magnetic transition at T 2 ~ 24 
K. In the polar phase below T 2 , a coplanar helicoidal 
magnetic structure is stabilized with propagation vector 
k rn = (0.4876,0.2813,0.2029) in a general direction in 
reciprocal space. Strikingly, the electrical polarization 
in Cri3Nb 2 08 (with a magnitude of 17.8 //Cm" 2 ) is al- 
most exactly perpendicular to the plane of rotation of 
the spins, in clear contradiction with the predictions of 
the cycloidal multiferroics model. We conclude that the 
electrical polarization in Cu 3 Nb 2 08 must arise through 
coupling of the chiral component of the magnetic struc- 
ture with the crystal structure, rather than the ordinary 
polar component — a novel mechanism that is only al- 
lowed in a small number of paramagnetic point groups. 
This mechanism may also describe the magneto-electric 
coupling observed in other compounds reported to have 
P perpendicular to the spin rotation plane 

Polycrystalline rods of Cu 3 Nb 2 0g were prepared by 
the solid state route, with stoichiometric amounts of CuO 
and Nb 2 C>5 being treated at 950° C for 36 hours. Single 
crystal growth was then carried out in an optical floating- 
zone furnace (Crystal Systems Inc.), giving high quality 
samples with dimensions ~2x2xl mm 3 . 

The crystal structure of Cu3Nb 2 0s, not previously re- 
ported, was solved by single crystal X-ray diffraction us- 
ing an Oxford Diffraction, Super Nova instrument (A = 
0.7107 A). Further details and a set of refined struc- 
tural parameters are given in the Supplementary Mate- 
rial. Data collection and reduction was performed us- 
ing the Agilent CrystAlis Pro software packa ge, and the 
crystal structure was refined using FullProf ll|. A re- 
liability factor of i?(F)=2.94% was achieved when fit- 
ting 691 independent reflections that met the signifi- 
cance criteria, / > 2.0er(7). Cu 3 Nb 2 0s was found to 
be isostructural to triclinic Cu 3 NbTaOs [II], [l3j]. There 
are two distinct Cu 2+ sites in the unit cell, located at 
Wyckoff positions la (Cul, on an inversion center and 
in square-planar coordination) and 2i (Cu2, in general 
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FIG. 1. (Color online) (a) and (b): The crystal structure 
of CU3ND2O8 in the ac- and 6c-plane, respectively. Cu 2+ 
ions form saw-tooth chains parallel to the a-axis with steps 
of three edge sharing CuC>4 squares. The steps are linked via 
Cu2-0-Cu2 "risers". Nb 5+ atoms separate the chains along 
the 6-axis. (c): The magnetic structure of the Cu 2+ chains. 
The envelope of the spin rotation is shown by grey circles, (d) 
P (green arrow) perpendicular to the spin rotation plane. 



position and square-pyramidal coordination) of space 
group PI. The magnetic Cu 2+ sites (S=l/2) form saw- 
tooth chains parallel to the a-axis with 3-atom Cu2- 
Cul-Cu2 steps sharing the edges of their CUO4 squares, 
linked by Cu2-Cu2 "risers" (jumps), sharing the edges 
of the triangular faces of the pyramids (figures QJ, and 
Q})). The chains are then separated along the &-axis 
by a layer of non-magnetic Nb 5+ ions. The room- 
temperature lattice parameters of Cu3Nb20g were re- 
fined to be a=5.1829(5)A, 6=5.4857(7)A, c=6.0144(7)A, 
o=72.58(l)°, /3=83.421(9)°, 7=65.71(1)°. 

Figure [2^ shows the temperature dependence of the 
specific heat of Cu3Nb2 0s, measured using a Quantum 
Design, PPMS. Two anomalies are evident, indicating 
phase transitions at 26.5 K and 24.2 K. Magnetic suscep- 
tibility measurements, obtained using a Quantum Design 
MPMS, confirmed the magnetic nature of the phase tran- 
sitions (figure HJd). At higher temperatures, the magneti- 
zation is dominated by a broad feature centered at 40 K, 
characteristic of short range magnetic correlations. A 
clear anomaly was observed at 26.5 K, concomitant with 
the first anomaly in the specific heat. Figure[2j; shows the 
electric polarization measured in four directions; (0,1,0), 
(3,2,3), (5,2,-5) and (6,2,6). The ferroelectric transition, 
labeled T2, coincides with the 24.2 K anomaly in the 
specific heat. This behavior is similar to many of the cy- 
cloidal multiferroics [H, [H, [H , in which an electrical polar- 
ization develops at the lower transition. The polarization 
vector, P, was determined to be in the (1,3,2) direction 
with a magnitude of 17.8 /iCm" 2 at 10 K. The four polar- 
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FIG. 2. (Color online) (a) The specific heat of Cu 3 Nb 2 8 . 
(b) The magnetic susceptibility of a single crystal sample, 
measured in a field of 1000 Oe || (3,2,3). The inset shows 
the same measurement up to room temperature, (c) The 
electric polarization in three approximately orthogonal and 
the general (6,2,6) directions, determined through the inte- 
gration of a pyroelectric current, measured at a warming rate 
of 1 Kmin" 1 , having poled the sample with E — 2 kVcm -1 . 
The (6,2,6) data was measured down to 21 K and extrapo- 
lated below, (d) The temperature dependence of the neutron 
magnetic diffraction intensity of the fundamental reflection at 
d « 10.4 A. The data has been fitted with a power law. 



ization components were calculated at 20 K (colored dots 
in figure [2k) and found to be in excellent agreement with 
the data. It is important to note that the electric polar- 
ization was 100% switchable in electric field, as shown in 
figured for the (0,1,0) direction. 

A powder neutron diffraction experiment, performed 
on the WISH instrument Q, ISIS, UK, confirmed the 
existence of incommensurate magnetic ordering below 
Tn = 26.5 K. Comparison of diffraction data mea- 
sured at 1.6 K (figure [3]) and 30 K clearly showed a 
number of well-correlated magnetic diffraction peaks, 
which can be indexed using the propagation vector 
k m = (0.4876, 0.2813, 0.2029) (incommensurate in all 
three reciprocal space directions). Additional diffraction 
patterns were measured at a series of temperatures be- 
low 30 K. The magnetic propagation vector was found 
to be approximately constant throughout the antiferro- 
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FIG. 3. (Color online) Rietveld refinement against neutron 
powder diffraction data of CU3ND2O8 at 1.6 K. Top and bot- 
tom tick marks indicate the nuclear and magnetic peaks, re- 
spectively. The measured and calculated profiles are shown 
with dots and a continuous line, respectively. A difference 
curve (obs.— calc.) is shown at the bottom. Inset: tempera- 
ture dependence of the b* component of k m . 

magnetic phase, except for a pronounced change in the 
b* component between 24 K and 26 K (inset of figure 
13]), close to the second magnetic transition at T^. The 
strong magnetic diffraction peak at d » 10.4 A (the fun- 
damental reflection where the scattering vector equals 
k rn ) was integrated to give the temperature dependence 
shown in figure [5Ji, which is proportional to the square 
of the magnetic order parameter. Power-law fits to the 
temperature dependence of the magnetic diffraction in- 
tensity (continuous line in figure and electrical po- 
larization yield critical exponents of (3m — 0.25 ± 0.02 
and Pp — 0.354 ± 0.001 (compatible with pseudo-proper 
magneto-electric coupling [15j). 

An initial fit to the neutron data with a collinear mag- 
netic structure wrongly predicted Cu (S = 1/2) mo- 
ments to be > 1 fiB- An alternative scenario of a ro- 
tating spin order is therefore clearly favored. However, 
it was not possible to find a unique solution from the 
powder data alone. Consequently, a non-resonant mag- 
netic X-ray diffraction (NRMXD) experiment was per- 
formed at beam line 116, Diamond Light Source, UK. 
This technique exploits the dependence of the NRMXD 
cross-section on the direction of the magnetic moment 
with respect to the incident and scattered directions of 
the light and the X-ray polarization [lti ]. To minimize 
the fluorescent background, the incident x-ray beam was 
tuned to 7.835 keV, well below the copper and niobium 
K absorption edges. Three magnetic diffraction peaks 
were then located. At each reflection the sample was 
rotated about the scattering vector (azimuthal scans) 
whilst measuring the diffraction intensity in both the un- 
rotated {a — a') and rotated (a — %') polarization chan- 
nels, employing a graphite analyser crystal tuned to the 
(006) reflection. The NRMXD data are shown in figure 
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FIG. 4. (Color online) Non-resonant magnetic x-ray diffrac- 
tion intensity measured as a function of azimuth angle at 
the a) (-3-k a * ,-2-k b * ,-3-k c *), b) (-3-k a * ,-1-fct* ,-3-fe c * ) and c) (- 
4+fc a * ,-3+fc;,* ,-4+fc c « ) reflections. The line shape, dependent 
on magnetic moment direction, has been fitted (red line) to 
find the spin rotation plane. 

21 The magnetic moment directions were evaluated by 
fitting the azimuthal dependence of the NRMXD cross- 
section, which was found to be extremely sensitive to the 
spin rotation planes. An initial fit with independent rota- 
tion planes for the two sites (Cul and Cu2) converged to 
a single plane within the experimental error. The normal 
vector of the common spin rotation plane {6 and 0), and 
scale factors for the two sites (Fi and F2), then gave four 
free parameters to be refined. In spherical coordinates 
defined with an orthonormal basis (xyz) such that a \\ x 
and b is in the rry-plane, the normal vector of the rotation 
planes was found to be at 9 = 75.5(2)° and cj> = 54.9(2)°, 
which corresponds approximately to the (1,2,1) recipro- 
cal space direction. The scale factor ratio, equivalent to 
the ratio of the squares of the magnetic moments, was 
found to be F2/F1 = 1.00(3). The discrepancies between 
calculation and data at the extremes of the a — n' az- 
imuthal dependences are likely to be due to systematic 
errors inherent to the measurement. Regardless, this er- 
ror in the intensity scaling does not affect the evaluation 
of the spin rotation plane. It may, however, lead to an im- 
precise determination of Fi / F± , which indeed differs from 
that determined by the neutron diffraction measurement 
described in the following paragraph. 

The spin rotation plane, as determined from NRMXD, 
has been constrained in the final Rietveld refinement of 
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the neutron powder data, enabling us to determine the 
relative phases between the copper sites, the ellipticity 
of the spin rotation and the absolute value of the mag- 
netic moments. The result of the refinement is plotted 
in figure El and has reliability factors R p — 6.46% and 
Rwp = 5.52%. The copper spins were found to describe 
a circular modulation (no ellipticity), with a magnitude 
of 0.89(2) fi B and 0.69(1) /j, B , on sites Cul (la) and Cu2 
(2i), respectively. The phase of the two ions occupying 
the 2i sites (inequivalent in PI symmetry) were refined 
to be 1.03(7)7r rad and 1.05(7)7r rad relative to the ion at 
site la. The magnetic structure, described in full in the 
Supplementary Material, comprises an essentially ferro- 
magnetic coupling within the 3-atom steps of the Cu 2+ 
chains, and a predominantly antiferromagnetic coupling 
through the jumps, as shown in figure [TJ;. The propaga- 
tion vector describes a generic helicoidal structure, with 
a fast rotation of the spins in the b* and c* directions, 
and a slow rotation parallel to a*. The most important 
result of our study was obtained by combining the results 
of the magnetic structure determination with those of the 
pyroelectric current measurements, and is illustrated in 
fig. [T](d). The electrical polarization is nearly perpendic- 
ular to the plane of rotation of the spins (~ 14° to the 
vector normal to the spin rotation plane), in clear contra- 
diction with the predictions of the cycloidal multiferroics 
model. In the remainder, we will show that our result is 
still compatible with the inverse DM model, but calls for 
an altogether different interpretation. 

Energy minimization in cycloidal magnets leads to the 
creation of local DM vectors Dy , associated with pairs of 
magnetic atoms i and j separated by the bond vector r y . 
In general, we can decompose Dij in components paral- 
lel and perpendicular to as = P^ x r,j + cr^Tij, 
where Py is a polar vector and cr^- is a pseudoscalar (both 
quantities are time-reversal even). Py can be taken to 
be the local polarization, so that the macroscopic polar- 
ization is P = Y^ij Pij- Conversely, a = a ij i s a 
macroscopic chirality. For perfect cycloidal and helical 
structures, <Xy = and Py = 0, respectively (a cycloidal 
structure is not chiral). For a generic helicoidal struc- 
ture, both Oij and Py are present, but only the latter, 
which is always within the rotation plane of the spins, 
is considered in the cycloidal multiferroics model. The 
chiral term <x,j cannot in itself produce a polarization, 
but it can do so through coupling with the crystal struc- 
ture, provided that the paramagnetic phase supports a 
macroscopic axial vector A, so that P — ^oA (7 being a 
purely structural coupling constant). Note that the elec- 
trical polarization in this case is always parallel to A. 
Macroscopic axial vectors represent collective rotations 
of one part of the structure with respect to another, and 
are only allowed in certain crystal classes, which we shall 
refer to as ferroaxial classes. There are seven non-polar 
ferroaxial classes: 1, 2/m, 3, 4, 6, 4/m, 6/m. Cu3Nb 2 0s 
belongs to the 1 class, the only one allowing a structural 



axial vector, and therefore an induced polarization, in an 
arbitrary direction. Another previously unrecognized ex- 
ample of a ferroaxial multiferroic is RbFe(Mo04)2 [13], 
(space group P3, crystal class 3) which has a proper he- 
lical magnetic structure. In RbFe(Mo04)2, the electrical 
polarization is again perpendicular to the plane of rota- 
tion of the spins, and is constrained by symmetry along 
the 3-fold axis. A more intuitive way to describe this phe- 
nomenon is the following: it is well known that the direct 
DM effect produces a proper screw magnetic structure 
in structurally chiral MnSi [l7|. Conversely, a magnetic 
structure with a screw (helical) component will induce 
structural chirality by the inverse DM effect. If the para- 
magnetic group is ferroaxial (i.e., supports a macroscopic 
rotation described by an axial vector A), an electrical po- 
larization P will develop in the direction of A. The sign 
of P is determined both by the direction of A and by the 
sign of the magnetic chirality. We note that in addition 
to the inverse DM effect, other more general microscopic 
mechanisms [3, 0, may apply to CuaNbaOg, which 
are also compatible with ferroaxial coupling. 

In conclusion, we have shown that in Cu3Nb2 0s an 
electrical polarization perpendicular to the plane of rota- 
tion of the spins develops below T2 ~ 24 K. We interpret 
this phenomenon, which cannot be explained within the 
framework of the cycloidal multiferroics theory, as due 
to the coupling between a macroscopic structural rota- 
tion (which is allowed in the paramagnetic group) and 
the magnetically-induced structural chirality. 
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